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We construct K7t amplitudes that fulfill exact elastic unitarity, account for one loop Chiral Perturbation 
Theory contributions and include all 1 /Nc leading terms, with the only limitation of considering just the lowest- 
lying nonet of exchanged resonances. Within such scheme, the Nc dependence of a and p masses and widths 
is discussed. Robust conclusions are drawn in the case of the p resonance, confirming that it is a stable meson 
in the limit of a large number of QCD colors, Nc- 

Less definitive conclusions are reached in the scalar-isoscalar sector With the present quality of data, we 
cannot firmly conclude whether or not the Nc = 3 /o(600) resonance completely disappears at large Nc or it 
has a sub-dominant component in its structure, which would become dominant for a number of quark colors 
sufficiently large. 
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Cl|' I. INTRODUCTION 

^ ' Light y^ — 0+ scalar resonance properties are of great interest, since they might help to unravel details of QCD chiral symmetry 

^S^ , breaking and confinement. Despite many theoretical efforts, the current understanding of the microscopic structure of these 

resonances is still far from being complete. The difficulty is triggered by the fact that scalar mesons carry vacuum quantum 

,__! . numbers, and also because strong final-state interactions hide their underlying nature when they are produced. Among other 

J> resonances, the Hghtest one /o(600), cuiTently denoted as the a meson, is an essential ingredient of the nuclear force, as 

t — anticipated long ago |1]. Its contribution to the mid-range nuclear attraction provides saturation and binding in atomic nuclei. 

■^!t| ■ During many years, there has been some arbitrariness on the "effective" scalar meson mass and coupling constant to the nucleon, 

partly stimulated by lack of other sources of information. The existence of this broad low-lying state is by now out of question; 

its mass and width have accurately been extracted from data analysis incorporating a large body of theoretical and experimental 

C^^ constraints |2-4]. The debate on the nature of the a meson is nonetheless not completely over. Structures of the tetraquark or 

^^ glueball type have been proposed (see e.g. Ref. [5] for a recent review and references therein). It is remarkable that such an 

^~~^ accurately determined state is so poorly understood from the more fundamental point of view of the underlying QCD dynamics 

. . [ of Nc^ Nf quarks and A^^ — 1 gluons, where A^c = 3 is the number of color species and Nf is the number of flavors. To clarify 

K^ ■ the issue on the nature of the a meson, it has been suggested to follow the dependence on a variable number of colors Nc 7^ 3 

of its mass and width IQ] by assuming that hadronic properties scale similarly as if Nc was large. A prerequisite for this scaling 

approach to work is that at least all leading-Nc effects are taken into account. 

d The limit of an infinite number of quark colors keeping UsNc fixed (ttj is the strong coupling constant), turns out to be a 

very useful and simplifying starting point to understand many qualitative features of the strong interaction ||7|,[a]- While keeping 

essential properties of Quantum Chromodynamics (QCD), under the assumption of confinement, the large-A^c limit provides a 

weak coupling regime to perform quantitative QCD studies. In this work we are interested in describing the large-A/c scaling 

of nn scattering and the induced large-A^ behavior of two-pion resonances. At leading order in \/Nc, the meson-meson 

scattering amplitudes are given by sums of tree diagrams induced by the exchange of an infinite number of weakly interacting 

physical (stable) hadrons. Indeed, meson and glueball masses scale as &{Nc) whereas the widths do as &{\/Nc) and G[l/Nc), 

respectively. Crossing symmetry implies that this sum is the tree-level approximation to some local effective relativistically- 

invariant Lagrangian, which by assuming that confinement still holds at large Nc, can be re-written in terms of mesonic fields and 

hence complies with quark-hadron duaUty. Higher-order 1 /Nc corrections correspond to hadronic loop diagrams and effectively 

restore unitarity in the time-like region. 

Resonance Chiral Theory (R;i;T) [l^l includes the pseudo-Goldstone bosons and the resonances as dynamically active 
degrees of freedom of the theory. The low-energy limit of R;tT must comply with low-energy theorems based on Chiral Pertur- 
bation Theory (ChPT) 1 1 1'-tO'], and this property has been used to predict systematically the Low Energy Constants (LECs) of 
ChPT in terms of masses and couplings of the resonances, when integrating them out of the action, at the chiral orders ff{p^) IMI 
and ff{p^) [14]. The ChPT Lagrangian includes the octet of pseudo-Goldstone bosons, however, when extending ChPT, R;i;T in- 
corporates the resonances as active degrees of freedom that are included in nonets, since octet and singlet of a SU [Np = 3 ) group 
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FIG. 1: Leading large-A'c TtTt scattering diagrams made of contact (polynomial) and resonance-exchange (pole) terms in t, s and u channels. 

merge into a nonet for Nq — > oo. The ChPT Lagrangian is built using the spontaneously-broken chiral symmetry of massless 
QCD. The explicit symmetry breaking from non-zero quark masses and electromagnetic interactions is incorporated in exactly 
the same way as it happens in QCD. The nonets of resonances are added requiring the general properties and invariance under 
charge conjugation and parity, and the structure of the operators is determined by chiral symmetry. At first order in the XjNc 
expansion, terms with more than one trace and loops are suppressed. The first property permits to postpone some terms allowed 
by the symmetries as subleading. However, the theory determined by just symmetries does not share yet some of the known 
properties of QCD at high energies or accepted from hadronic Regge phenomenology^. Further constraints on the couplings 
arise by matching the interpolating resonance theory at intermediate energies with asymptotic QCD at the level of Green func- 
tions and/or form factors. The application of these properties determines a series of relations between the couplings of R;ifT, 
reducing the number of couplings and enhancing the predictability. We remind here that QCD, besides current quark masses, 
has only one dimensionful parameter, Aqcd- Of course, there are infinitely many such short-distance constraints and hence a 
sufficiently large number of states may eventually be needed to avoid contradictory results llTsll . 

A further and less trivial question is related to whether or not purely contact terms should be regarded as independent of 
exchange terms in the large-A^c framework. Within %% scattering, this corresponds to distinguish between contact \% vertices 
and those where a resonance field propagates between 2n states (see Fig.[T]). While the large-A^c expansion is expected to provide 
a better approximation to data in the space-like region where unitarity does not play an active role, the approach is not specifically 
related to a given energy range. This means that if the tower of all infinitely many states is included, the question on the duality 
between contact and exchange terms is pertinent. However, if the exchanged resonance spectrum is truncated above a given 
energy, contact terms must necessarily arise to encode the explicitly disregarded high-energy contributions while still complying 
to the short-distance constraints 1 1011 . The explicit values of the contact LECs depend on the functional parametrization adopted 
for the resonance fields. 

In the Single-Resonance- Approximation (SRA) scheme, each infinite resonance sum is just approximated by the contribu- 
tion from the lowest-lying meson multiplet with the given quantum numbers. This is meaningful at low energies where the 
contributions from higher-mass states are suppressed by their corresponding propagators. The SRA corresponds to work with 
a low-energy Effective Field Theory (EFT) below the scale of the second resonance multiplets. In this work, we will use the 
SRA of RjT to re-analyze how the C7 and p properties depend on Nq, mainly based on the study of %% scattering. An early 
investigation was proposed in Refs. |'16',47'] keeping the leading 1 J'Nq contributions but omitting the ChPT chiral logarithms. We 
impose leading-A^ short-distance constraints, as discussed in detail in Refs. 1 10, 18]. This turns out to be of capital importance 
because it leads to a clear distinction between leading and subleading Nq contributions to the n% amplitudes. It also allows 
for a meaningful extension of the framework to the Nc > 3 world. Besides, enforcing the short-distance constraints reduces 
the number of free parameters to only the subtraction constants, needed to restore exact elastic unitarity, and the masses of the 
exchanged resonances. In addition, we check that a direct analysis of ;r;r scattering at leading order in l/Nc and using high- 
energy constraints, based on forward dispersion relations and Regge phenomenology, generates relations between the resonance 
properties compatible with those akeady obtained in Ill0l[l8ll by looking at other processes. 

We will first construct nn amplitudes that fulfill exact elastic unitarity, account for one-loop ChPT contributions, and include 
all 1 /Nc leading terms in the SRA. Next, we will look for poles in the appropriate unphysical sheets of the amplitudes, and 
discuss their properties when Nc deviates from its physical value," to learn details on their nature. In this manner, we improve 
on previous analyses lla.l20l-l22t] where leading 1/Nc terms, beyond a certain order in the chiral expansion, were neglected. 

On the other hand, in the strict chiral and large-A^ limits, the pseudo-scalar singlet 771 and the tt are degenerate ll23l - l25ll . 
The interplay between ChPT and large Nc has been addressed in Ref. |26]. U{3) meson-meson scattering has been treated in 
II27I1 with only contact 6'{p^) chiral interactions. In Ref. [,28,1 . in addition to the contact iff{p^), the leading-At scalar resonance 
exchange has been included, but leaving aside the vector meson exchange. This latter mechanism not only contributes to the 



' Regge behavior, while not directly deduced from QCD, works rather well and it has played a decisive role in the benchmark and unprecedented accurate 

determination of the a meson mass and width |2-4], by extending the high energy region above \/s > 1.4 GeV. 
^ We will focus here on the properties of the poles in the complex s plane. The Nc behavior of the Breit-Wigner resonance parameters has been discussed, in a 

model independent manner, at length in Ref. Lift! , and we refer the reader to this latter work for further details. 



P-waves, but also to other S-wave channels via the left-cut contribution. Recently, a full-fledged one-loop unitarized couple- 
channel t/(3) calculation, including both scalar and vector resonance exchanges, has been undertaken in Ref. II29|] . Besides 
achieving an excellent description of phase shifts until center-of-mass energies of around 1 .4 GeV, this reference also analyzes 
the Nc behavior of the amplitudes. Among other results, it is explicitly shown there that when Nc increases, the mass of the 
lowest eigenstate of the rji-rjn mass matrix decreases, reaching values of around twice the pion mass in the Nc — 30 region. 
Thus, a natural question arises here, namely, do the T7-T7' degrees of freedom play a relevant role to determine the Nc trajectory 
of the (7 and p resonances? This is addressed also in Ref. f29], from where one might infer that this is certainly not the case. 
Those degrees of freedom turn out to be much more relevant in the study of the Nc dependence of masses and widths of higher 
resonances, as for instance the /o(980). We benefit here from this observation, and we will neglect tj-tj' effects in what follows. 

We should also point out that some aspects of the A^ ^ 3 extension undertaken in Ref. |29] deserve discussion, and we 
believe they can be improved along the lines followed in this paper In particular, the leading 1/Nc contributions were not 
properly considered in Ref. |29]. As a consequence, we cannot firmly conclude a scenario where the a moves far away in the 
complex plane for large Nc, as obtained in 1.29.1 . We will give some more details below Eq. ( |56] |. Note that crossing symmetry, 
the possible absence of an exotic isotensor (/ = 2) state and the well-established fact that the p width decreases with Nc imply 
the existence of a narrow scalar-isoscalar resonance, in the large-A^c limit, with a mass comparable to that of the p meson [3Q0. 

The paper is organized as follows. In Section |II] we introduce our notation in order to make the paper self-contained. We 
also review some important features of nn forward dispersion relations and crossing symmetry, which yield to two sum rules, 
keeping an eye on the large-A^c expansion. In Section|III]we discuss the nn scattering amplitude in the SRA of RjT and analyze 
high-energy conditions which yield to the short-distance constraints. In Section [TVl the 6'[p'^) ChPT contributions to the kk 
scattering amplitude and its matching to the leading \/Nc piece are discussed. Actually, we manage to write an amplitude which 
is correct to &{p^) and includes, within the SRA, the leading ff{l/Nc) terms to all orders in the chiral expansion. However, it 
still needs to be unitarized before being confronted to scattering data. After unitarization, we fix our parameters in Section [Vl 
by fitting the (/,7) — (0,0), (1,1) and (2,0) phase shifts in the real A^c — 3 world. This complies with the determination of the 
R;i;T Lagrangian parameters at the leading-A^ approximation. Only then do we allow ourselves to analyze the Nc > 3 situation 
in Section[VT]and the emerging picture for nn resonances as the number of colors is varied. Finally, in Section lVlIl we draw the 
main conclusions of our work. 

11. GENERAL FEATURES OF ;r;r SCATTERING 



A comprehensive presentation of nn scattering can be seen at the textbook level 113111 and more recently in 113211 . We summarize 
here the relevant formulae to fix our notation and to provide a proper perspective of our subsequent analysis merging large Nc, 
ChPT and unitarity considerations. 

A. Kinematics 

The ndpi) + %(p2) -^ ^c(p'i) + ^d{P2) scattering amplitude is written as 

Tab-cd ^A{s,t,u)5ab8cd+A{fj^j'^)^acAd+Mu,t,s)5ad5bc, (1) 

with the standard choice of Mandelstam variables s = [pi +P2)~, t = {p\ — p[)^ and u = (pi — pj)^, where A {s,t,u) is the 
n^n^ — > Tr'^Tr" amplitude, which is the only independent amplitude thanks to isospin, crossing and Bose-Einstein symmetries. 
That is, if Tj^{s,t,u) is the isospin combination with total isospin / (in the ^-channel), one has, 

T,^=o{s,t,u) = -{3A{s,t,u)+A{t,s,u)+A{u,t,s)}, (2) 

Ti,=i{s,t,u) = -{A{t,s,u)-A{u,t,s)} , (3) 

T,^=2{s,t,u) = -{A{t,s,u)+A{u,t,s)}. (4) 



For the normalization we will use here the conventions from II33I1 . The partial-wave decomposition in the .^-channel becomes 

Ti{s,t,u) ^f^{2J+l)Tu{s)Pj{cosO), (5) 



where Tij{s) is the projection of the nn elastic scattering amphtude with given total isospin / and angular momentum 7: 
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being nijc = 139.57 MeV the pion mass and Pj the Legendre polynomials. The in-elasticity riij{s) = 1 for s < 16m^ and 
riij{s) < 1 for s > 16m^. Besides, 5/y are the phase shifts and the Mandelstam variables t and u depend on s and on 0, the 
scattering angle in the center-of-mass frame (cm.). The optical theorem reads 






where the partial-wave total cross section is defined as 

(27+1) 
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This value is bound by 



au{s)<%%{2J+l)^^-^. 
s ~ 4mi 



(8) 



(9) 
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The contribution of a resonance state, with spin J and isospin /, to the partial cross section in the narrow-width hmit and assuming 
m^j ^ 4m^ reads. 



mij 



(11) 



Of course, one may think that such a limit does not apply to a broad state as the a. Within a Breit-Wigner model, the finite 
width correction effectively corresponds to a reduction F ^ F [1 — Y/{nm)\, which even for the extreme case F = m yields to a 
moderate 30% correction. 

Using the recent GKPRY parameterizations of Ref. [34] for the partial S, P, D and F waves, cross sections are presented in 
Fig.|2]up to ^/s < 1.42 GeV. As we see, the SO, P and DO waves play an outstanding role featuring the appearance of the /o(600), 
Pi (770), /o(980) and /2(1270) resonances. Below the KK production threshold, s^^ = Am\, only SO and P are essential. 



B. Crossing and forward dispersion relations 



In the (crossed) /-channel the amplitudes read 

^7>,=o(i,f,M)^ 

r/,=i(i,f,M) 

,7/,=2(«,f,M), 
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Jh=2{s,t,u) 



(12) 



where 7, and 7^ are the corresponding isospins in the t- and ^-channels respectively. 

The rigorous Froissart bound from axiomatic field theory requires that in the forward direction (t = 0) these amplitudes, up to 
logarithmic corrections, are asymptotically polynomially bound by a single power of s. Actually, in terms of the crossing-odd 
variable v = {s-u)/2, the amplitude 7}, (v,f) = Ti,{v + 2ml-t/2,t,-v + 2ml-t/2) satisfies T},{-v,t) = (-)''7),(v,f) -^ v"'', 
with rtji < 1 when v ^- oo. This means that a forward (t — 0) once-subtracted dispersion relation is fulfilled, by considering a 
closed contour excluding the cuts —°° < V < — 2m| and 2ot^ < V < +°° (see Fig. |3]l, where T{v + /0+) — T{v — /0+) — 
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FIG. 2: (Color online) The total partial-wave cross sections Oij{s) (in GeV ), defined as a function of the cm energy variable yG (in GeV) 
for the SO (Solid, Red), DO (Solid, Black), P (Solid, Blue), F (Dashed, Blue), S2 (Solid, Brown) and D2 (Dashed, Brown) partial waves below 
^/s = 1.42 GeV, from the GKPRY parameterization with set UFD of parameters given in |,34t]. We also draw the unitarity bounds (dashed 
lines) for SO (red), P (blue) and D2 (black) waves, using t]oo = '?ii = 1 and T702 = 0-88 respectively. 
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FIG. 3: The contour in the complex crossing-odd v plane used to derive FDRs for the /-channel isospin kk scattering amplitude. For 
/ = V = .s — 2m^, so that v = 2m^ corresponds to the ;r;r threshold. The principal value of the forward scattering amplitude corresponds 
to 7/^(v + iO+). Finite Energy Sum Rules arise for the Regge- subtracted amplitude for the finite circle |v| < V/,. Low-lying resonances 
/o(600) = (7, pi (770), /o(980) and /2(1275) are also marked by crosses in the Second Riemann Sheet, across the cut. 



2ilmT{v + /0+) = 2/Imr(v). The subtraction constant can be fixed at low energies, and more specifically at v = or equiva- 
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FIG. 4: The /, = 0, 1 , 2 total cross sections, as a function of tiie s variable, including S, P, D and F waves below \/s = 1 .42 GeV (solid lines) 
and the corresponding Regge behavior starting at ^/s = 1.42 GeV (dashed lines) used in the solutions of Ref. 1341 . /< = (Red), /< = 1 (Blue) 
and /, = 2 (Black). 



lently s = 2m^, i.e. below threshold. Thus, at f = 0, we get the Forward Dispersion Relations (FDRs): 



r,,=o(v,0) = 7),=o(0,0) 
7>,=i(v,0) 



2v2 /■" c/v'Imr/,=o(v',0) 



= lim 

Vo^O 
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v'^ - v^ 



Vo 



fl'v'lmr/,=i(v',0) 



7},=2(v,0) = 7),=2(0,0) + ^ 



n Jlmi V'2 V'2- 

c/v'Im7),=2(v',0) 



The absorptive part of the amplitude can be written as an optical theorem in the v variable: 



Im7},(v,0) = -Jv2-4m4(J/,(v), 



(13) 



(14) 



where one goes from /, to /^ with the same s-t crossing matrix as for the amplitudes, see Eq. (fT2] i. The FDRs converge, provided 
|Imr/, ( V, 0) I < v" with a <2, for large values of v. Results for (7/, are presented in Fig. 21 up to ^/s = 1 .42 GeV, using the UFD 
parameterization of Ref. |34] for the partial S, P, D and F waves 

These general constraints are obeyed for any Quantum Field Theory and by QCD in particular. On the other hand, for f 7^ 0, 
and sufficiently large v, say v > V/,, one has the phenomenologically successful Regge behavior given by^ 



Im7},(v,0^-27r2^j3/,(r) 



«;,(0 



(15) 



where ^indicates summation over several Regge trajectories, which for small f are given by «/,(?) = «/, (0)+fa/ +. . . The slope 
parameter, a'j^ (0)^1/ {2m^ ) = 0.9 GeV"^ is nearly universal and the leading trajectory intercepts are Oo (0) ^ 1 , ai (0) ^1/2 
and 0(2(0) ^ (for a review within a modern Ttn context see e.g. Refs. 135113611 and references therein). We also show in Fig.|4] 



Our amplitude and that used {F{v,t)) in Ref [34J are related by T{v,t) = -27[^F{v,t) 



the Regge behavior used in the construction of the partial waves via FDRs, Roy and GKPY equations.'* As we see, there is some 
mismatch between the partial waves and the Regge behavior. This is somewhat expected, since Regge behavior provides an 
average of the oscillating resonance contribution in the high-energy region. The region ^/s < 1.42 GeV is well described by S, 
P, D and F waves, since the exchanged p meson in the f -channel corresponds to longest range 1 /nip and hence /max ^ \fslmp . 

C. Finite Energy Sum Rules 

Separating the Regge tail in the dispersive integral, corresponding to the integration region v' > V/,, we have for /, — 0,2, in 
the limit v -^ +°° and assuming a/, (0) < 2, 
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and similarly for It — 1, 
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(17) 



(18) 



with a = o:/,=i(0) — 1. Note that the last term in both Eqs. (1% and (fTSl) is a constant subleading contribution. We will denote 
this constant term as — /?/, . Requiring 



^^^(v,0) - ^-'L-^^-"''^'^^^'(iy'"\ ^^=0,2, 



Ti,=iiv,0) ^ 2K^y Pi,=m ^-.a.n ( y 



(19) 

(20) 



in the v — > +°° limit, we get the Finite Energy Sum Rules (FESRs) 
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(21) 



** For our purposes we need the following Regge tails |34], valid for .v > .v^ = (1.42 GeV)^: 

Im7},=o(v,0) = -2K'^[hp(vlvo)+bp{vlvofi'\., 
Im7),=i(v.0) = -2Tl^bi(v/voY\ 
lm7},=2(v,0) = -2;r2fo2(v/vo)°2. 

The numerical values of the parameters are given in section 8 of the appendix A of Ref. 1341 . 



(16) 



where 



\VqJ ap, \Vq 

AKb, 1 /v,^"' 
fll-1 V/, Vvo 

R,.2 =An'-l('A"\ (22) 

For ^ = V V/, + 2m2 = 1 .42 GeV, we get the values 7},=o(0) = -121.8+64.3 (P)+28.2 (P') = -29.3, 7};=j(0) = -105.2- 
19.4 (p) GeV-2 ^ _i24.6 GeV^^ ^^^^j 7)^^2(0) = -12.25 +7;/,=2- Indeed, in the window 0^ G (1, 1.42) GeV, we observe a 
smooth V/j dependence of both the integral and the Rj^ contributions.^ Moreover, the total sum itself remains fairly independent 
of V/,, as well. Being more quantitative the right-hand side of Eqs. ^T\\ changes at the level of 10% and 4% for the /, = and 
1 cases, respectively, when ^/vj, varies in the interval (1,1.42) GeV. For the case /, = 1, R'j ^i amounts, at energies as small 
as 1 GeV, to around 25% of the total, and decreases with increasing V/,, as deduced from its Regge behavior. For the case of 
/, = 2, the integral contribution shows a more pronounced dependence on V/,, and it changes sign at around half of the (1,1.42) 
GeV interval. However, this contribution in size is just at maximum around one third of that of /?/,=2, which is almost constant, 
because of the smallness of the Regge intercept, 02- Furthermore, the sign of Ri,=2 depends on the sign of 02 and actually the 
FESR becomes ambiguous for 02 — 0, since the amplitude would show a logarithmic growth ^ log(— v^/V/^) instead of a Regge 
behavior. For a diverging/converging amplitude we have a positive/negative contribution to the sum rule. These ambiguous signs 
appear in the CFD and UFD solutions of Ref. 1,34] , from where one gets Ri,=2 = —30(30) and ^/,=2 — 20(40), respectively, both 
results compatible with zero. 

D. Adler and a Sum rules 

If/?^^j and/?/,=2 are neglected in Eqs. (l2Tl i. as suggested by the numerical values obtained from the fits performed in Ref . 03411 . 
and Tj ^^{Qi) and Ti^jjO) are approximated by the lowest-order result in ChPT, see Eq. (|52] | below, the so-called Adler and a 
sum rules (see e.g. Olll for a discussion based on current algebra and PCAC) are obtained. 



1 _ 2 n dv i 4»4 
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3(^0(V) + -CT1(V)--C72(V) 
3Cyo(v)--C7l(v) + -C72(v) 



(23) 
(24) 



with /;/i: ^ 93 MeV, the pion decay constant. Using the GKPRY parameterizations Il34ll . one finds that the contribution of the 
region v > 1.42 GeV, for which Regge behavior is assumed, is quite small. This is because (T/,(v) — > 471^)3/, (O)v"'''-'''^^ In 
addition, we would like to point out 

• /; = 1 : Following the discussion of the previous section, neglecting the contribution to the sum rule arising from the region 
\/v G (1, 1.4) GeV might induce variations of order 20%, which we expect to be of the same order as those stemming 
from ff{p'^) ChPT terms neglected in the left-hand side of the sum rule. 

The Adler sum rule is satisfied to 5% (see also IBTT ) in the partial-wave plus Regge representation. 

• /, = 2: This sum rule converges if a/,=2(0) < and is the a sum rule derived in [31] on the basis of PCAC. The o sum 
rule is already approximately satisfied for an upper limit of the integration around ^ < 1 .26 GeV. The contribution above 
this upper limit tends to cancel and shows an oscillating behavior as a function of the upper limit of the integration. This 
is closely linked with having a value of ^2 compatible with zero and a not well-defined sign for the ratio ^2/02 appearing 
in R/,=2- For an upper limit of y/s < 1.26, the truncated Adler sum rule also provides a rather reasonable value of 89 
MeV for f„. This result is certainly more reasonable if one bears in mind that the left-hand side of the sum rule has been 
computed at lowest order in ChPT only. This observation suggests a kind of super-convergent dispersion relation which 
will be important to set up our model below. 



' We use the UFD parameterization of Ref. 
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Given the above discussion, one finds reasons to saturate the sum rules with the lowest-lying resonances below 1 GeV. In 
addition, we already mentioned that the I, —Q sum rule in Eq. ( ITTT i is saturated also, with great accuracy, at 1 GeV (integral and 
R contributions tend to cancel above 1 GeV). All this, gives support to the scenarios that will be presented below, in which we 
will take into account only the lowest-lying resonance degrees of freedom. Yet, we will make use of the Adler and cr sum rules, 
saturated at energies of about 1 GeV, to find out short-distance constraints that will make more predictive the Rx^ approach 
adopted here. 

III. LARGE-Nc ASPECTS OF nn SCATTERING 

In this section we discuss some aspects of the large-At limit which will provide some useful guidance in our analysis of nn 
scattering. 

A. 7tK elastic scattering amplitude in the SRA 

From the lowest-order R;)fT Lagrangian l^llfll, describing the couplings of the lowest-lying V{1 ) , A(l++), 5(0++) and 
P{0 ^) resonance nonet multiplets to the pions, we find* 



A'"''^{s,t,u) 



s Gy (t{s — u) u{s — t) 



+ 



fn /k \t-ml u- ml 

2 [c,i{s~2ml)+2mlc,„] 



2 



3/4 s - ml^ 

4 [cd{s-2ml)+2mlc,„] 

tA 2 

Ml ml 

In the large-A^c limit, |q/| — |cj|/a/3 and |c,„| = |c,„|/-\/3. We have specified for clarity the contributions from non-degenerate 
singlet. Si, and iso-singlet octet, ^8, fields. Quite generally, msg ~ms^ — 0{\/Nc) and mixing effects have been analyzed in 
Refs. mUltl.^ Glueball mixing within Rx'^ has been discussed in Ref. [40]. 

Taking m^^ = m^, = m^, |q/| = |cj|/\/3 and |c„,| = |c„,|/a/3, we reproduce the expressions in Ref. IHlll . In principle, the 
couplings appearing in the scattering amplitude can be determined by analyzing the decay processes p ^2n and S ^^2n with 
S ~ ■\/2/35i + Ss/V^, corresponding to a (uu + dd)/\/2 flavor composition in the qq picture, which yield in the chiral limit 
(see also Ref. jil).^ 



s — 



'' ^ (26) 



lenf; 



•4 






* Here we use the antisymmetric field formulation where the A-71 mixing is absent. Note that the axial resonance does not contribute to the elastic mz scattering 
amplitude. After proper incorporation of short-distance constraints, the Proca formulation using gy = Gy / fit yields the same amplitude |10]. (/(S) nonet 
resonance fields are generically parametrized as R = ^ Y^=i ^i^i + ^^o with A; the standard Gell-Mann matrices. 

' Ref. 13911 finds sizeable mixing effects between the singlet and octet scalar-isoscalar mesons and identifies two possible phenomenologically acceptable 
scenarios for the tree-level mass eigenstates: a) Ml =1.35 GeV andM^ = 1.47 GeV, or b) Ml =0.985 GeV andM// = 1.74 GeV. The light solution was less 
preferential at it would correspond to a case where /o(980) would not couple to pions. Therefore in several studies the heavy solution has been adopted. 

Clmi _i , r, ^ 7 / „2 



' The full expressions are Fj- = , ^ pj 



Crf + (Cm-Cd)-? 



^"'* ^^ = li^Pv' *'''=■■'= P^ = Vl"4m2/„ 
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The residues of the scalar-isoscalar and the vector-isovector poles in the partial-wave amplitudes are 



gs 



■2 



J K 



Note that the (large-A'c) relations g^ — 167irsms/3 and gy = 16;rrymy hold in the SRA, in the chiral limit. The SRA amplitude 
in Eq. dZST l contains too many parameters to be analyzed in full detail. In the next subsections, we will discuss a sensible way of 
reducing the number of couplings and masses. 

After projecting onto partial waves we get for the {I, J) — (0,0), (1, 1) and (2,0) channels the following asymptotic behavior 
at large values of s: 



^SRA/„^ -i 2q + 3Gy - f^ 5 1 a ^ a ^ 

71 S+ ... , Aoo — 1 , All — -, ^1-20 — - T 

Ik 6 2 



Tjf\s)^Xu ''^•"r/ '^ s+..., Aoo-1, All--, A2o = --. (29) 



The above behavior implies that subtractions would be necessary to make convergent a dispersion relation.^ We will add 
subtraction constants after unitarization. 

B. Short-distance constraints 

The short-distance constraints encompass R;i;T with proper high-energy behavior 11811. In general they produce a set of 
conditions which, for a limited set of resonances and in particular in the SRA, reduce the number of independent couplings. 
These conditions may be over-determined yielding at times to mutually inconsistent values, a problem which can be side-stepped 
by introducing more resonances. In the present case we think it of interest to pursue such an analysis within ;r;r scattering. 

In the chiral limit, the f-channel amplitudes, see Eq. (fT2b . in the forward direction have the following asymptotic behavior in 
terms of the crossing-odd variable v = {s — u)/2: 

r-A(v,o) - 2«±%i±^M + ^(v-), 

In 



^SRA,„n^ _ 6cj-3fl+4gl + 3Gl 
3/^ 



T,f^tiy,o) = -. "■'^^r/^'-^ v+^iv-'), 



r-,^v,0) = 2fM±M^Z^MZ! + ^(v-), (30) 

where we have included momentarily the tensor meson coupling gr and mass nir, to be discussed below in more detail. This 
limit is compatible with the Froissart bound, a specific merit of the antisymmetric tensor formulation [KT].''' 

As already mentioned, it makes sense to divide the nn scattering amplitudes in Eq. ( fT3] l into three pieces: i) the low-energy 
part which takes the form of subtraction constants and is fixed by chiral symmetry, ii) an intermediate-energy part, dominated by 
resonance exchange, and iii) the high-energy remainder which we expect to be responsible for the Regge behavior. Therefore, if 
we impose a behavior for the resonance contribution no worse than suggested by Regge theory, we obtain the constraints 

fn = 2cj + ^gl + Gl, (31) 

= 6clmj + 4-glml-3Glml. (32) 



' Requiring that the leading term proportional to A/j vanishes would yield the relation 2q + 3Gj, = /^, advocated in Ref. |42]. This would imply Gv < fjc/V^, 
giving a 30% too small p ^ 271 decay width. This relation has also been found as a necessary high-energy constraint in a NLO R;|;T analysis of the vector 
form factor, incorporating subleading l/Nc corrections |43]. In the Q = Hmit, it gives Gv = /k/V^, which was also found in the study of one-meson 
radiative tau decays carried out in 1 44j] . Notice however that imposing this relation is not enough to make subtractions unnecessary, because the partial waves 
would still grow at large values of .v as T/j (.v) ~ (mj, Gy //^ ) log {s/ifiy ) . 

Physical results are actually independent of the field representation. The naive exchange of Proca fields does not satisfy the Froissart bound, but after suitable 
polynomial subtractions to comply with the short-distance constraints, one ends up with the same amplitude llOil - Fields remain a useful framework to 
incorporate symmetries, see also the discussion in Ref. |45j. 
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These constraints correspond to require that the v and the v'' coefficients of 7)^5f and 7)^^^, respectively, vanish in the large-v 
regime. The second condition is less robust than the first one, attending to what we have discussed on Regge phenomenology 
above. Indeed, Eqs. dsTl i and ( |32] ) can be also obtained from the Adler and a sum rules, in the chiral limit, using the narrow 
resonance approximation of Eq. ( fTTI ) to estimate the cross sections that appear in the right-hand sides of the sum rules. Note, 
that since we have not considered any exotic isotensor resonance, we are approximating (72 = 0. 

In the absence of tensor couplings, gj- = 0, these constraints imply c^ = sm^f^/V^, Gv ^cos^/jr and my /rris = tan0, where 
is a mixing angle. The KSFR relation (Gv ~ fji/^/2) requires ^ = 7r/4 and hence nis = my, as well as 2£j = \/2Gv = fn, 
implying F^ = 9Tv /2. These constraints have also been found in the algebraic chiral-symmetry approach i46l - l48ll and can be 
re-written, in terms of the decay widths, as 



r^ 32Kfl ^ 9 Ty 32nfl 
nis 3m| 2 my 3my 



1 = ^_^ + --L_4£, (33) 



^ ^-1^, <34) 

nis 2 my 

which yield the value of the scalar mass and width to be 

ms = 660 MeV , Tg = 570 MeV , (35) 

when phenomenological values for the mass and width of the p meson are used. These numbers are quite sensitive to details. 
For instance if the KSFR set of parameters is used, one gets instead (taking as input the p meson mass) 

ms = my/ = 775 MeV , Fj = ^Fy = 805 MeV . (36) 

Of course, given the fact that the scalar turns out to be a broad resonance, it is unclear what these estimates should be compared 
to, since generally a resonance is characterized by the complex pole and the complex residue of the scattering amplitude. The 
benchmark calculation of the pole on the second Riemann sheet of the kk scattering amplitude [2, 3], when written as Sa = 
m^ — im<jF(j, yields m^ — 347(17) MeV and Fa = 690(48) MeV. On the other hand, the connection between the Breit-Wigner 
(BW) resonance parameters, defined as 5(mg^) = k/2 and Fbw — l/('MBw5'(m|\y)), and the pole resonance parameters has 
been discussed on the light of their A^ behavior in a model-independent fashion 1 19], suggesting that the large shift in the mass 
is ff{\/N}-) and can be computed, yielding an acceptable extrapolation of mew ^ 700 MeV. The recent TTTT-scattering analysis 
of Ref. (ImI leads to the Breit-Wigner values [m^, Bw,r<j, bw] = [841(5) MeV, 820(20) MeV].'' The quoted errors above also 
account for the existing differences when UFD and CFD parameterizations are used. This yields a ratio F^ bw/'W(j. bw ^ 
5.0(1) Fp^ Bw/'«p, BW, which suggests a 10% accuracy of large Nc in the SRA, supporting as well the identification of the 
large-At parameters with the BW ones. 

C. Higher-energy resonances 

We have so far been limited to states below the KK threshold, ^/i < 1 GeV On the other hand, Regge behavior works for 
y/s > 1 .4 GeV So, it is interesting to see the modifications induced by other resonances which may decay into 2n, in the mass 
range 1 GeV < m^ < 1.4 GeV, namely /o(980), /ij (1170), ^i (1235), /2(1275), /o(1370) and pi (1450). 

After implementin g th e appropriate short-distance constraints, via the Froissart bound, the inclusion of a 2++ tensor yields a 
resonance amplitude ll49ll 

AT{s,t,u) = -^^ ^ '- 11^ 2 -^ (37) 

/^ mj--s f2 mj. 

where the coupling is determined from the decay into ;r;r in a relative D-wave yielding 

2 3 



" The value of the pole is s/J^ = 445(8) — 1297(7) MeV, In agreement with Ref. |2]. On the other hand, the model-independent large-A^ -based extrapolation 
from the resonance pole mass to the BW pole mass 1 19] yields m^, bw = 670(20) MeV, when the UFD parameterization of Ref. |34|| is used. 
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The previous amplitude yields the following contributions to the 6'{p'^) ChPT couplings 11491] : L^ = Lj = and L^ = 
gj/{3mj) ^ 0.16 X 10^^ (see also Ref. pO']). The axial 1+^ mesons, such as /ii(1170) and /7i(1235), give just a purely 
polynomial contribution to the nn scattering amplitude (without .^-channel propagator poles) which cannot satisfy the Froissart 
bound, yielding to no contribution at all to the LECs. We remind that the exchange of 7 > 1 resonances in the f-channel naively 
violates the Froissart bound, a situation which has been the standard motivation to rely on high-energy Regge behavior as a way 
of introducing suitable cancellations. 

If the tensor meson /2(1275) is considered, we might include p' = pi (1450) and /o(980) as well, where the decay widths into 
Tin are taken to be T{f2 — > TtTz) = 150 MeV, r(/o — > Ttn) — 80 MeV and T{p' -> TtTt) — 250 MeV (note the lai-ge inaccuracies). 
This yields the extended sum rules (we remain in the chiral limit) 

^ ^ y Ts 32Kfl ^ y9Yv l,2nfl , y g^r 327r/^ 
■^ nis 3ot| y 2 my 3my j mj 3m^ 

~ nig r^ mj y 2 my 

Using PDG values fsl^, the higher resonances /o (980), /2( 1275) and pi (1450) produce corrections of the order of (0.06, -0.01) 
for the r.h.s of the first and second sum rules respectively. This shows a trend to cancellation which supports that higher states 
not only play a minor role at low energies but also in the region of interest below 1 GeV, to leading order in Nc- Therefore we 
will carry our analysis below with just scalar 0++ and vector 1 states.'^ 

D. Other short-distance constraints 

Alternatively to the previous analysis, one may derive short-distance constraints from other processes involving two- and 
three-point functions Ill0l[l4ll . Imposing the short-distance properties of the underlying QCD dynamics, within the SRA, one 
gets ladSJ: 

V2Gv = 2c,/ = 2c,„ = 2\/3 cj = 2V3 c„, = 2\/2ii„, = f„ . (40) 

These constraints are obtained from a variety of processes, some of them involving electroweak probes. It is remarkable that 
they turn out to be totally compatible with those deduced here by looking to Ttn phenomenology at high energies, for a vanishing 
tensor-meson contribution (gj- — 0). Neglecting the /2(1270) and higher-mass resonance contributions leads to a realistic and 
simplified scenario for the purpose of the present work. On the other hand, from the ;r;r scattering amplitude, we have derived 
the additional restriction ms = my ■ Although we will explore the effects of this constraint in one of the fits that will be presented 
below, we should mention here that it relies on the assumption of an asymptotic behavior for 7/,=2(v) more convergent than that 
of a constant (v"). This is not a totally robust result, though it is certainly plausible, given the accuracy of current nn analyses 
at high energies |34]. 

For my = 775 MeV, the conditions in Eq. ( |40l i lead to 

r. = ^p. = 750MeV, (41) 

3 

^y = 5?72Pv = 150MeV. (42) 

Besides, by requiring the two-point correlation functions of two scalar or two pseudoscalar currents to be equal at high energies, 
up to coiTections of the order Usf'^/t^, one finds [18] mp^ — V2ms. This relation involves a small correction, of the order 5%, 
which we neglect, together with the tiny effects from light quark masses. We will use Eq. ( l40b to fix our parameters below, but 
allowing ms to be different from my . 



'^ Of course, it is intriguing to analyze tlie role of the (fi ,j'''~) = (2+,0^^) exotic state, named X(1420) in the PDG |51] (and found also in the SU(6) study of 
Ref. ^2\), which decays into 2;r in a S wave with estimated mass Mx = 1420(20) MeV and width Tx = 160(20) MeV. The contributions of such a state to 
the Adler and fj sum rules, Eqs. )23l24t , in the narrow-resonance approximation are — SOTrPx/Sm^ ~ —4.6 GeV"~ and IdnTxI^nix ~ 1.9, respectively. As 
compared to the individual contributions from other mesons, they seem too small to provide a clear signal. 
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IV. ONE-LOOP CHPT CORRECTIONS AND UNITARIZATION OF THE TtTt ELASTIC SCATTERING AMPLITUDE 

A. TtK scattering at one loop in CliPT 



At ff{p ) in ChPT, the TtK elastic scattering amplitudes can be written in the form Ill2ll : 

A'^^^'^{s,t,u) = A^^^'^{s,t,u)+A'^^^'^{s,t,u), (43) 

A?"(.,r,«) = 4r^, (44) 

In 

AjhPT(,^f_„) ^ __i_|(2[i +[2_ Z),2 + (/-2_ £)(f _«)2+4(3;;_2/"i - ^)mls - {3h + UU - 8h - y>4 

+ 3 (,s^ - ml)J{s) + (t{t-u)- 2m\t + Am\u - 2m\) J{t) + {u{u-t)~ 2m\u + Am\t - 2»4) /(«) | . (45) 

The lowest-order ampUtude A2''^^(s,/,m) is identical to the first term in Eq. ( [25l l (pion contribution) and only depends on the 
pion mass and weak decay constant. The ff{p'^) correction involves four SU(2) renormalization-scale-independent LECs: U 
(/ = 1,2,3,4). In addition, A^'^^^ (,?,?,«) includes one-loop chiral corrections, which are suppressed by one power of l/Nc', they 
are parameterized through the loop function 

-p{s)-l- 



J{s) =2 + p{s) log 



PW + 1 



(46) 



For the sake of completeness, we recall here, the relation between the LECs /,■ and the most common SU(3) parameters L'-{pL) 
OB: 

32;r2 ^ ^ 

I, = — -/;-(M)-log(m2/Ai'), (47) 

where 

7i==T' 72 = T' 73==-^, 74 = 2 (48) 

and 



1 2 

-3' ''-r 


73== 


1 

^2 


l[ = 4L[+2Li- 


"24' 





V, = -8^-4L^ + 16L|: + 8L^-^, 

Vv 

II = 8Li + 4L^-^, (49) 

with 327r2vA:,7) = l+\og{m^^/^^),m-q = 4otjij'/3 and m^: ^ 468 MeV the kaon mass in the limit m„ = m^/ = 0. The renormalized 
coupling constants /[(/x) and L\{pL) depend logarithmically on the dimensional regularization scale pi: 

Lj(M2)=L;(Mi) + Y^log(g), (50) 

where 

3 3 1 3 11 5 

El = — , F, = — , F^ = . F4 = - , F5 = - , Eg = , Eg = — . (51) 

^ 32 16 ' * 8 8 144 ' 48 

The corresponding f-channel forward scattering amplitudes at v == are given by: 



7>,=o(0) 



m 



2 »4 (-72/1 -48/2 + 45/3 + 36/4 + 397: -101) 



2/2 5767r2/4 

1 (-48/4 + 9;r-1 0)m2 ^ ^^ g^ 
/2 + 384;r2/4 



n-M - -i+^^^^^^S7r^+^(/)' 



r,.2(0) ^ -f + (-48/2 + 18/3^72^ + 21.+ 10K ^ ^^^.^ ^^^^ 

Ik JibK Jjf 

The ^{p'^) corrections are at the 10 — 20% level of the lowest-order ones in the Adler and a sum rules. 
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I ChPT-improved SRA amplitudes and large-Nc counting rules. 

In the limit of a large number of colors QISl,^^'^ in Eq. dZSl ) scales like l/A^c since the pion weak decay constant behaves 
like ff{y/Nc). Furthermore, A^^^ provides the leading-A^c prediction for the actual nn scattering amplitude, with the only 
limitation of considering just the lowest-lying nonet of exchanged resonances |9]. This latter approximation is justified as long 
as s, t and u are kept far from the second resonance region. 

The lightest resonances have an important impact on the low-energy dynamics of the pseudoscalar bosons. Below the reso- 
nance mass scale, the singularity associated with the pole of a resonance propagator is replaced by the corresponding momentum 
expansion; therefore, the exchange of virtual resonances generates derivative Goldstone couplings proportional to powers of 
l/m| . At lowest order in derivatives, this gives the large-A^ predictions for the ^{p'^) ChPT couplings |9]. At ff{p'^) the 
amplitude A^*^^ in Eq. dZSl) reduces to 



Af^(.,f,M) = 



m 



2 



s Gy (t{s — u) u{s — t) 



fl fn I ml ml 

2 [cd{s-2ml)+2mlc,„] 

4 [cd{s-2ml)+2mlc,„] 
ft mj^ 

dj ml 
f4 m^ ' 



2 



+ ^4, (53) 



which constitutes the leading l/Nc approximation to A*^''^^. The polynomial form can be re-written as 

°'^^?r f /Of SRA I rSRAN^ , /ajSRK . o,SRA /i r SRA orSRAN„^2\ /^hx 

~ ^^1(2^4 +^5 J'^+l'^i-e +2^8 ^4L4 -2L5 )m^\ . (54) 



From Eqs. ( l53T l and (|40] |. one trivially finds 



r SRA _ r SRA 



In J SRA _ -'In , Jk 



Of SKA rSKA J n j \ 

8m^ ttniy 8mj 

Lf^^O, if'^-^, if^^O, ^r^^lfl, (55) 

in full agreement with Ref. lUSl] .'^ 

Let us pay now some attention to the Nc dependence of the one-loop ChPT amplitude. Note that from Eq. (l45T l. the logarithmic 
contribution to A^''^^ scales as 1 jN^, while the polynomial piece behaves as 1 JNq in the A'c ^ 1 limit. This is because, as shown 
in the relations ( fSSl l. the LECs L, behave as ff{Nc), with the exceptions of L2 — 2L\, L4 and Lg that scale as ff{N^) [13]. The 
renormalization scale dependence of the LECs provides further subleading contributions in the I/Nq counting. Unfortunately, 
the measured values of the L,- couplings cannot be phenomenologically split into their large-A^c leading and subleading parts. In 
general, one has the scale-dependent relation 

mii)=AiNc + Bi{}i), (56) 

where A,- is scale independent. Note that only the Nc = 3 combination is experimentally accessible. However a meaningful 
extension of the chiral amplitudes to an arbitrary number of colors requires some knowledge of the coefficients A, and B,. This 



If the relation ms = my is furtlier assumed, one gets purely geometrical ratios for the non-vanishing LECs: 



TrSRA_,SRA_ ^ ,SRA _ ' rSRA _ ^ sRA _ fit 
2Li -L. ---L3 --L5 --L, -^ 



In chiral quark models with meson dominance built-in one also has ms = my = \j2AfNcnf ~ 781 MeV, for / = 88 MeV |53], so that the above relations 
hold with L2=A'c/192;r2 = 1.5 X lO^^ 
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difficulty is precisely what Af^^ helps to overcome, and thus the SRA predictions of Eq. dSST l can be used to read off the A, 
coefficients.'^ 

In the recent work of Ref. |29], the SRA parameters Gv, cj, c,„, Cd,Cm, my, nisg and m^j are fitted to data. Afterwards and 
to extrapolate to Nc > 3, the resonance masses are kept constant, while all the couplings are scaled by a y/Nc/3 factor. This 
does not take into account that singlet and octet scalar resonances become degenerate in the large-A^c limit. On the other 
hand, if one fits the resonance couplings to data, the fitted values do not necessarily match their leading-A^c values in Eq. ( |40] l 
and might incorporate some significant l/Nc subleading contributions, which later on however, are scaled as if they were 
leading in the Nc counting. For instance, in Ref. |29] a value of around 15 MeV is found for Cd, which is around a factor 
three smaller than that of fj[/2 quoted in Eq. ( l40l i. A proper extension of this parameter when Nc deviates from 3 should 
be 2cj = A*^^ X y/Nc/3 + (30 MeV — /^i*^ ), instead of that assumed in [29]. For other parameters, there exist also large 
deviations between the fitted values found in |29] and the leading-A^c estimates given in Eq. ( |40l ). These differences seem to 
indicate 1 /Nc corrections to the resonance parameters which are much larger than expected. One might wonder the underlying 
origin of these disturbing large deviations. The fitting strategy certainly might play some role on this; for instance, the choice of 
the upper energy limit or the choice of the unitarization procedure. This latter issue has some relevance which we will address 
now. In the next section, we will discuss our unitarization procedure, which is rather similar to that used in Ref. ||29I1 . There 
appear independent subtraction constants for each of the {I, J) — (0,0), (1,1) and (2,0) sectors Il33ll54]. However, in Ref. 112911 . 
all three subtraction constants were forced to be equal. From the discussion in 133115411 . this is somehow an arbitrary choice. 
The lack of flexibility of data fits incorporating such constraint might influence the actual values determined for the resonance 
couplings and their estimated scaling with Nc- 

Let us come back to our scheme. The amplitude A^*^^ in Eq. (|54] | contains neither pion loop terms, nor the 1 /Nc subleading 
contributions to the polynomial piece of the one-loop ChPT amplitude. Thus, and to better describe the experimental phase shifts 
for Nc = 3, we propose to use 



A^'^+^"^'(i,f,M) = A''^is,t,u)+ A^"'^Hi,f,M)-Af^(i,f,M) 



.ChPT/ . „\ ;.SRA/ 



(57) 



e{i/Nc) 



In this way, by construction, we recover the one-loop ChPT results, while at the same time, all terms in the ampli- 
tude that scale like l/Nc (leading) are also included, within the SRA. Note, that in the l/Nc counting, the correction 
I^ChPT (^ f [<) _ A4^^(i, f , m) } is incomplete, since it does not account for all existing subleading 1 /N^ contributions to A{s,t,u). 
A complete l/N^ calculation would require quantum corrections stemming also from the low-lying resonances 11431 155145911 . 

C. Unitarized amplitudes 

The ChPT loops incorporate the unitarity field theory constraints in a perturbative manner, order by order in the chiral expan- 
sion. Though subleading in the l/Nc counting, their effect appears to be crucial for a correct understanding of 5-wave nn phase 
shifts. Furthermore, resonances show up as poles of the unitarized amplitudes in unphysical sheets, which positions provide 
their masses and widths. Thus, to discuss the nature of the /o(600) resonance, it is crucial to restore unitarity. Any unitarization 
method re-sums a perturbative expansion of the scattering amplitude in such a way that two-body elastic unitarity, 

lmTj]'{s) = ^, s>4ml, (58) 

is implemented exactly. To restore unitarity, we will make use here of a once-subtracted dispersive representation of Tij^{s) 
(see for instance Sect. 6 of Ref. ||^). Let be .^jj'^^'^^^'^ (s) and [T2]ij{s), the TZTt ampHtudes, in the {I, J) sector, deduced from 
^SffA+ChPT^^ J j^-j and Aj'^^^ («,?,«), respectively. We define an unitarized amplitude as li33il 

7)7 1(.) = -Cjj~Io{s)+Vr/{s), (59) 

Vij{s) - ^i'j'^+'''''^{s)^[T2]ij{s)iIois)+Cjj)[T2]ij{s), (60) 

where Vij stands for the two-particle irreducible amplitude, and the subtraction constant and the loop function are given by 

CiJ = -7)7' (4^2 ) + Vrj\4ml), Io{s) ^J^i^- -/»] • (61) 



''' Actually, the SU(2) scale-independent LECs defined by Eq. j471 display the Nc separation in a scale-independent fashion, since Ij = ajNc + b where b is 
common to all coefficients and stems from pion loops. This allows to build differences, _// — Ij = [ai — aj)Nc, which can be used to extract the leading-A^ 
contributions from data up to a constant. This is illustrated in quark model calculations 15311 . 
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Note that i) the constants C/y determine the scattering length/volume in each sector and become undetermined free parameters, 
and ii) with the election of Vjj^{s), and considering A^'^'^(i,f, m) ~A^^^{s,t,u) as ^{p^), we recover from T[j{s) the ChPT series 
up to one loop, in all {I, J) sectors. 

Analytical expressions for the ChPT amplitudes projected onto the different {I, J) sectors can be found in the appendix B of 
Ref. il. 

A final remark concerning the l/Nc counting rules is in order here. The subtraction constants Cij must scale as 

Qj^mi^), (62) 

because of their definition as the difference between the inverses of the full and the two-particle-irreducible amplitudes at 
threshold for each {I, J) sector.'^ Thus, the amplitude Tjj{s) reduces to the 1/Nc leading part of 3^ij ^ (s) in the limit of a 
large number of colors. 

V. Nc = 3 PHASE SHIFTS AND p AND fo(600) MESON PROPERTIES. 

In this section, we fit the Qj parameters to phase-shift data and show results for the poles found in the second Riemann sheet 
(SRS) of the amplitudes. The SRS of the T matrix is determined by the definition of the loop function Io{s) in the SRS. We 
follow here Ref. [60] and use Eq. (A13) of this latter reference to compute 7o(.s') in the SRS. '^ 

Masses and widths of the dynamically generated resonances are determined from the positions of the poles, sr, in the SRS of 
the corresponding scattering amplitudes in the complex s plane, namely sr = M| — i MrTr. For narrow resonances {Tr <C Mr), 
y/sR ^ Mr — irfl/2 constitutes a good approximation. '^ 

The coupling constants of each resonance to the pion pair are obtained from the residues at the pole, by matching the ampli- 
tudes to the expression 

for values of s close to the pole. The couplings gR are complex in general, and represent independent information from the pole 
Sr. In the narrow-resonance approximation, the extrapolation of Eq. ( |63] ) to the real axis takes a Breit-Wigner form to comply 
with unitarity and hence g^ = \6nmRT r / p {m\) . 

The first issue is to select the set of data points to be fitted. We will consider the scalar-isoscalar, vector-isovector and 
scalar-isotensor elastic nn phase shifts, with a total of 107 data points, as follows. 



• I = J = sector: As our main input, we will use the Roy-equations results from Refs. yl |61|, |62|] in the energy range 
v^ < 750 MeV. We take this upper cm. energy cut to keep negligible coupled-channel ^/T effects, which give rise to the 
/o (980) resonance. We have considered the phase-shift determinations of Ref. 161!] and that of Eq. (4.8) of Ref. |62] . For 
each value of y/s, we have used as central value the average of both results, while the absolute difference between them 
is taken as the error for the x^ fit- From threshold to the upper cut of 750 MeV, we have moved up in steps of 10 MeV, 
which amounts to a total of 48 phase shifts to be fitted. 

• I = J=1(I~ 2, J = 0) sector: We fit to the phase shifts compiled in Refs. ll63ll64ll ( ll65ll66ll ) and consider an upper cut 
of V^ < 910 MeV (1190 MeV), which comprises a total of 38 (21) phase shifts. 

An important point has been the selection of the fitting interval. Clearly we should restrict the range to the elastic region and, in 
particular, below the opening of the first inelastic KK channel. At lowest order, this effect corresponds to the two-step process 
nn — > KK — > nn, which in the sub-threshold region yields a real contribution to the amplitude and is l/N^ and i^^^-suppressed. 
The pure elastic re-scattering is just 1 /A^^-suppressed. On the other hand, corrections to the SRA due to heavier states with 
mass M are of ff{l/{NcM^j). Therefore we expect that, by restricting to low energies, inelastic effects can be safely included 
in l/Nc corrections to the LECs. In other words, we may allow M ^ 2Mk ^ 1 GeV without much trouble. 



'^ The full and the two-particle-irreducible amplitudes differ by terms which always contains at least one i-channel pion loop, which is subleading in the l/Nc 
power counting. Thus, the difference T^' — V^' scales as ff{N^), as one readily deduces by noting that both TV and {T — V) scale as ff{\/N^) in the N^ S> 1 
limit. 

'* Note that the function Jo[s), defined in Eq. (A8) of f6ul, reduces to h{s) for equal masses. 

'^ It has become customary to quote for the broad fj the number ^Js^ instead of s^. Here, we will also quote ^/s^, but we will refrain from identifying Rey's^ 
with the mass and — llmy^s^ with the width of the resonance. For the p meson, the narrow-resonance approximation works within two standard deviations. 
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The second issue is to design the fit procedure. For three out of the five fits examined here, we will fix 

mv = 0.77 GeV, ms = 1 GeV. (64) 

For the fit B.3, we will force my = ms, as suggested from our discussion of Eqs. dsTl i and (l32l ). and fit the common value to data. 
Finally, in the fit B.1-2 we will fix my to 0.77 GeV, while ms is considered as a free parameter We will always fit Coo, Cn and 
C20 to data; besides those parameters, my and ms — mp/^/l, ^■^^+ChPT depends still on Lj 2 3 4 5 g ^(m), once the relations in 
Eq. (|40l ) are implemented. 

We have considered two well-differentiated scenarios: 

• A: Wejust take the SU(3)Gasser-Leutwyler parameters LJ'(/i), at a certain scale ji, from other phenomenological analyses. 
In this type of fits, only the C/j parameters are fitted to data. 

• B: The contributions of the low-lying vector, axial-vector, scalar and pseudoscalar resonances to the L,-, and therefore 
to the effective chiral Lagrangian at order ff{p'^), were given in Eq. dSSl) . and thus, the renormalized coupling constants 
L[(/x) can be written as a sum |19|] 

mii)^L^^+mii) (65) 

of the resonance contributions, L?^^ and a remainder U-ill). The choice of the renormalization scale /i is arbitrary. 
However, it is rather obvious that one can only expect the resonances to dominate the i[(/i) when /i is not too far away 
from the resonance region. Therefore, it is common to adopt n = mp as a reasonable choice. However, one might take as 
a best fit parameter one scale, jJ.^^, for which it occurs a complete resonance saturation of all the LECs L- , this is to say, 

^■(m"^)=0. (66) 

In other words, at this privileged scale fi'^^, there is no other contribution in addition to the meson resonances. In this type 
of fits, besides the Qj parameters, this scale should be fitted to data. 

We have also considered a scenario where the complete resonance saturation of the LECs LJ occurs at two different 
scales, Hy^ for L\ ^ ^nd \i^^ for ^4 5 g g, depending whether the LEC is dominated by the vector or the scalar resonance 
contribution. Note that, L3 is renormalization-scale invariant. 

Let us start discussing results obtained from five different fits to the phase-shift data. Best fit parameters and pole properties are 
compiled in Table H] while predicted phase shifts and SRS amplitudes are depicted in Figs.|5]|6] and|2l In the first fit (A), we 
take the SU(3) Gasser-Leutwyler parameters L\{\i = mp) from the G{j>^) Ke4 fit compiled in Table 2 of Ref. Il67ll .^^ A best fit 
to data using, instead, the central values of the main ff{p^) fit of Ref. 16711 leads to a X^ jdof more than twice larger than that 
of fit A. This is not entirely surprising, since we are only considering here one-loop chiral logarithms, besides those required to 
restore exact unitarity, in the %n amplitudes. In the four other fits (B.l, B.2, B.1-2 and B.3), we assume that complete resonance 
saturation occurs, the corresponding scale (or scales) where it holds is fitted to data. In the fit B.1-2, we fix /ly'^ to the value 
obtained in fit B.l and we fit \i^^ and ms to data. Besides in the fit B.3, we force ms — my in the SRA amplitude, and fit this 
common mass to the data. We explore this possibility because it is suggested from our previous discussion on short-distance 
conditions in Eqs. OTT l and (l32b . The ms = my constraint was deduced in Ref. 1,221 as well, when the one-loop SU(2) ChPT 
amplitude, unitarized with the Inverse Amplitude Method (lAM), was required to be consistent with the SRA. 

Since in the SRA amplitudes we have explicitly incorporated one vector an one scalar poles, we expect at least one pole in 
each of the I — J ~ and I — J ~ I sectors. Because of the re-summation in Eq. i5% . the pole positions will change with 
respect to those of the bare ones (s = my and s = m|) and the resonances will acquire a width that accounts for their two-pion 
decay. In addition, as we will see, some other poles are generated as well in the SRS of the scattering amplitudes. The five fits 
have reasonable values of x^/dof, though in the 7 = 7=1 channel B-type fits lead to a better agreement with data than fit A. 
The major improvement occurs at the higher end of the fitted region, and it is due to the the tail of a second resonance located 
at around 1.4 GeV (see Fig. [T) that it is generated in the B schemes. The relation g^ = ml/ (6/^), deduced from the KSFR 

prediction Fp = mi/{967:f^), is satisfied within 3% accuracy. This is because the chiral logarithms are almost negligible in the 
p-meson channel, and indeed at order ffip^) they cancel out for the SU(2) massless pion theory ||22I1 . 

Besides, it is also interesting to compare the leading-A^ values of the residues, as displayed by Eq. ( |28] |. with the correspond- 
ing ones after unitarization, given in Table|T] Using the input values one gets \gy \ ~ 2600 MeV, in excellent agreement with the 



'* In units of 10 ', ^i 2 3 4 5 6 8(^ = "ip) = 0.46, 1.49, —3.18, 0., 1.46, 0. and 1.08, respectively. For simplicity in the analysis, we have ignored the eiTors on 



^l,2,3,4,5.6,8VH — '"p; 
these parameters, since they do not affect the main conclusions of this worlc. 
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FIG. 5: Theoretical predictions for the phase shifts obtained from fits A, B.l, B.2, B.1-2 and B.3. Fitted data from Refs. 163116411 (/ = \,} = 1) 
and 163 . [6ql (/ = 2,7 = 0) are also displayed. In the isoscalar-scalar (/ = 0,7 = 0) channel, the used average of the results of Refs. I6ll[6al is 
shown, as well. 



values for \gp \ quoted in the table that range in the interval 2400 — 2500 MeV. For the scalar resonance the rather stable values 
for \ga I ^ 3000 — 3 100 MeV can only be reproduced by Eq. ( |28] ) if m^ — my (Fit B.3). Indeed, for the B.3 fitted nig — my mass, 
we find \gs\ ^ 3000 MeV. For completeness, let us mention that the residue of the a resonance has also been deduced from 
the (7 ^ 77 decay [68], yielding ga = 3204(28) + i 1588(14) MeV, i.e. \ga\ = 3580(30) MeV. The result from Roy equations 
yields \ga\ ^ 3300(300) MeV |69]. Using the UFD and CFD parameterizations from lH yields gs = 3735(61) +1874(3) MeV 
and ga = 3742(60) +1874(6) MeV or, equivalently, l^^l = 3836(85) MeV and \ga\ = 3843(84) MeV respectively, while in the 
most recent analysis based in GKPY equations, a value of |^fy| = 3590(120) MeV is quoted [4]. 

Moreover, the scalar-isoscalar phase shifts are significantly better described when the complete resonance saturation of the LJ 
occurs at two different scales fitted to data (fit B.2). In this latter case, mass and width of the /o (600) or a resonance compare also 
well with the results of Caprini et al. 121 -'^ The properties of the <y resonance are strongly influenced by chiral logarithms [22!]. 
In this scalar-isoscalar channel, we will keep track of a second resonance, which cannot be identified with the /o(600) and that 
we will label as scl. There exist two types of scenarios: i) for fits A, B.l and B.3, this second pole appears well above mg (1 
GeV for the first two fits and 0.738 GeV for the last one) and it is relatively narrow, and ii) for fits B.2 and B.1-2, it is placed 
below nig (I GeV and 1.295 GeV, respectively) and it is quite wide (F > 600 MeV). In the case of the fit B.2, the effects of 



' Bear in mind that in the / = / = () sector we do not fit directly the Roy equation results of Ref. 1 61], used in the work of Caprini et al, but rather we take an 
average of these results with those obtained by Yndurain and collaborators in Ref. |62]. The recent re-analysis 1341 provides errors fully compatible with this 
assumption. 
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A 


B.l 


B.2 


B.1-2 


B.3 


Coo 


-0.0210 (5) 


-0.0218 (9) 


-0.0278(18) 


-0.0283 (18) 


-0.0144 (5) 


Cu 


-0.02054(11) 


-0.01996(11) 


-0.01979 (12) 


-0.01968 (12) 


-0.0108 (9) 


C02 


-0.0594 (19) 


-0.0588 (20) 


-0.0621 (17) 


-0.0593 (19) 


-0.0606(18) 


Ai^^ (MeV) 


770 (fixed) 


693 (26) 


474 (16) 


693 (fixed B.l ) 


520 (21) 


IX§'^ (MeV) 


770 (fixed) 


Mf 


1550(180) 


1190(130) 


Mf 


ms (MeV) 


1000 (fixed) 


1000 (fixed) 


1000 (fixed) 


1295 (40) 


738 (3) 


my (MeV) 


770 (fixed) 


770 (fixed) 


770 (fixed) 


770 (fixed) 


ms 


X^dof 


3.5 


2.4 


2.0 


2.0 


2.2 


contrb. ;t2 


25/323/21 


26/203/22 


5/178/20 


17/171/20 


12/193/20 


ri2,ri3,n4,ri5 


0, 0,-,- 


0.06,0.13,0.80,- 


0.04, 0.05, 0.43, 0.20 


0.09, -0.01, 0.47, -0.05 


0.20,0.08,0.51,-0.20 


r23,r24,'-25 


0,-,- 


0.01,0.07,- 


-0.02,-0.12,0.21 


-0.01, 0.29, 0.27 


-0.02, -0.20, -0.99 


I'34,r35 


-,- 


0.16,- 


0.14,-0.10 


-0.03, -0.03 


0.15,0.02 


'■45 


- 


- 


-0.59 


0.60 


0.22 


mp (MeV) 


749.1 (4) 


752.6 (5) 


754.0 (5) 


754.8 (6) 


754.6 (7) 


Tp (MeV) 


144.5 (3) 


150.4 (4) 


152.0 (4) 


153.1 (5) 


149.5 (4) 


l^lp (MeV) 


2404 (4) 


2490 (5) 


2515(5) 


2534 (7) 


2504 (6) 


V^ (MeV) 


(451±2,-;-234±l) 


(453±4,-;238±4) 


(423±6,-!267±3) 


(442 ±4,-/248 ±2) 


(446 ±4,-/246 ±4) 


IslcT (MeV) 


3005 (21) 


3080 (70) 


3070 (120) 


3100(50) 


3070 (90) 


m.d (MeV) 


1340 (40) 


1600+»2oo 


772 (6) 


1020l^^'o 


1030 (50) 


r.d (MeV) 


1,7+22 
11/ „ 


300+^0" 


5801120 


1070+24" 


200+™ 


kl.d (MeV) 


2800 (300) 


4700+»,„o 


2980(100) 


2940 (120) 


3300+5"" 



TABLE I: Best fit parameters and pole properties (statistical uncertainties on these latter quantities define 68% confidence-level regions, 
induced by the corresponding Gaussian correlated errors of the different fit parameters). The three Cjj parameters are always fitted to data, 
and in addition /Xy = ^s, ^v and ^5, fi^ and ra^, and fly = /I5 and my = nig are also adjusted in the case of fits B.l, B.2, B.1-2 and B.3, 
respectively. In the row labeled as contrb., the contributions to the X of the different [I = J = 0)/(/ = J = l)/(/ = 2,7 = 0) sectors are 
displayed. Besides, r/j are Gaussian correlation coefficients between parameters ; and j. Note that the dispersive data analyses based in 
Roy [2] and GKPY ^equations predict for 0^ = (441+^*', -/272+^2) MeV and (457+[4, -i-279+^' ) MeV respectively, while the Review 
of Particle PropertiesBUl quotes mp = 775.49 ± 0.34 MeV and Fp = 149. 1 ± 0.8 MeV. 



this resonance (scl) on the phase shifts, in the higher end (600-750 MeV) of the fitting range, are appreciable and considerably 
improve the achieved description (see Figs. |5] and |6]l. A different question is whether or not such a wide state, below 1 GeV, does 
have a correspondence with any physical state or it is just an artifact of the fitting procedure. We should note that a state with 
these features has not been reported neither in the Roy equation analysis of Ref. |69], nor in the most recent work based in the 
GKPY equations of Ref. |4]. As we shall see in the next subsection, the Nc ^ 1 behaviour of the ff— pole obtained from the B.2 
fit is radically different to that inferred from the A, B.l and B.3 schemes. That is the reason why we have proposed the fit B.1-2, 
with the intention of testing to what extent the dependence on Nc of the cr- resonance properties is determined by the existence 
of this possible artifact^". In the fit B.1-2, where the scale jj.y^ is fixed to the result of fit B.l and the value of ms is adjusted to 
the data, the second scalar resonance shows up above 1 GeV, and as we will discuss below, it leads to qualitatively the same A^c 
dependence of the /o(600) mass and width as the fit B.2 does. 



VI. RESULTS FOR Nc > 3 



We extrapolate the amplitudes to Nc > 3, by means of the Nc dependence 

2 



, SRA+ChPT 



{s,t,u) 



Nc>i Nc 



A'^%s,t,u)l^_^+^^^ (AChi^(.,f,«)-Af^.,f,«)) 



Nc=3 



(67) 



^^ The minimization procedure, involving also the election of the parameters which are fitted to data, is not completely unique, and there exist obviously different 
local minima. Some of them might not have a proper physical interpretation. Bearing this in mind, we can not discard the existence of artifacts. 
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FIG. 6: Modulus of the SRS Tqq{s) amplitude, as a function of complex s (s = M^ — 'iMF) taken in the fourth quadrant, for fits A, B.l, 
B.2,B.l-2 and B.3. In all plots the physical scattering line (F = 0) is also depicted. 



and the scaling law of Eq. ( |62] |. Nc > 3 results are depicted in Figs.[8ll9land[T0l 

In what the p meson properties concerns, we observe that for all five fits examined here (Fig.©, both mass and width behave 
as expected from a qq picture. Thus rup, that did not deviate at Nc = 3 much from my, quickly approaches to my, while Fp 
decreases like 1 /Nc, as the number of colors increases. This is not by any means a new result, and in the past other authors have 
already reached, within an unitarized ChPT scheme, this conclusion OalZllj. Indeed, in a previous work [22], the same result 
was obtained starting from the one-loop SU(2) ChPT amplitude for massless pions and using, as in Refs. ||6|,|2l|], the lAM to 
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: M^ — IMF) taken in the fourth quadrant, for fits A, B.l, B.2, 



restore elastic unitarity. However, results here are more robust, because in previous works the leading l/Nc terms appearing 
beyond 6(ji^) [6, 22] or ff{p^) ll2lll were simply ignored. Note that the constraint niv = nis + ff{l/Nc), deduced in Ref. ]22] 
when the one-loop unitarized nn amplitude was required to be consistent with the SRA, does no longer necessarily hold here, 
and we could still have both parameters to be independent (fits A, B.l, B.2 and B.1-2). Though this could be because we keep 
here all l/Nc terms at all chiral orders, it might also happen that the above constraint was just an artifact of the lAM used in 
02211 to unitarize the amplitudes. Nevertheless, we should note that fit B.3, where the constraint my — nis is enforced leads to 
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phenomenologically acceptable results as well. 

Let us move on, and discuss now the scalar-isoscalar channel. For sufficiently large A^, and since unitarization corrections 
are subleading, we should end up with just the unique resonance (of mass nis) included in the SRA irreducible amplitude, while 
the effects of the second resonance must disappear. We see in Figs. |9] and [TO] that this is effectively the case. However, now it is 
difficult to draw robust conclusions and the resonance that survives depends on the fit procedure. 

In the case of the fits A, B.l and B.3, we see that the resonance identified as the /o(600) for Nc = 3 becomes the SRA pole. 
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with mass ms when Nc is sufficiently large. There is a first transition region for values of At close to 3, where both the mass 
and width increase with Nc, but above At = 6 — 8 the resonance width starts decreasing like 1 /Nc, and the mass approaches to 
the limiting mg value. The behavior showed in the two upper panels of Fig. |9]is almost identical to that of the right upper panel 
of Fig. 1 in the two-loop analysis of |21]. The authors of this latter reference conclude that in the case of the a resonance, 
there exists a mixing with a qq sub-dominant component, arising as loop diagrams become more suppressed at large Nc- The 
nature of the a resonance in the real world (Nc — 3) would be totally different to that of the p meson, being it mostly governed 
by chiral logarithms stemming from unitarity and crossing symmetry 112211 . justifying the widely accepted nature of the ff as a 
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FIG. 10: Same as Fig.jS] but for the scl pole. Now, poles appear in the third quadrant for fits A, B.l and B.3. 
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dynamically-generated meson. However, within this scenario, for sufficiently large Nc > 10, the structure of both (a and p) 
resonances is similar As pointed out recently in Ref. 1 30], this solves the seeming paradox of how a distinctive nature for the p 
and a at A'c = 3 is reconciled with semi-local duality at larger values of Nc- Semi-local duality requires the contribution of these 
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elastic cross section to cancel^' "on average", since this process is purely isospin 2 in the f-channel, 
and there are no isotensor resonances at low energies. 

In what respects to the second resonance found in the scalar-isoscalar channel, we see in Fig.[TO]that for fits A, B.l and B.3, 
it follows a totally different pattern with increasing A^^. Indeed, it becomes wider and wider and from one value of A^c on, the 
pole Ssci turns out to be located in the third quadrant, though y/s^d still lies in the fourth quadrant, with Re i < (where the 
path integral for the resonance field would not be well defined) and its effects on the scattering disappear. We would like to 
mention that in Ref. [29] poles in the SRS are being searched in the variable y/s. That could be inappropriate as we pointed 
out in Ref. |22], and we reiterate here. This is because, as mentioned above, there are situations where y^ lies on the fourth 
quadrant of the SRS, but however sr has passed to the third quadrant and thus its meaning becomes unclear. This phenomenon, 
which can only happen for broad resonances, was also illustrated in Fig. 1 of 122(1 . and it is precisely what happens for the a 
meson case in Ref. |29] (see the Fig. 10 in that reference). The conclusion of Ref. 02911 that the a meson moves far away in the 
complex plane for large Nc overlooks the fact that it does so in the third quadrant of the complex plane. 

The recent work of Ref. 1 70] might contradict the findings for the a resonance deduced from fits A, B.l and B.3, which in turn 
seem compatible with those obtained in the two-loop analysis of Ref. ll2lll . In OTOll . instead, emerges a picture more consistent 
with that outlined in the one-loop analysis of Ref. [6]. Within the unitarized quark model proposed by Tornqvist |71], the authors 
of 1701 find that the whole low-energy scalar spectrum below 2 GeV, except for a possible glueball/o(1710), could be described 
in one consistent picture, with the bare "qq seeds" dressed by the hadron loops. In this model, the a resonance is generated as a 
pole of the S matrix and has no correspondence with any of the bare qq seeds included in the scheme. Indeed, the <7 resonance 
runs away from the real axis on the complex s plane when Nc increases. However, following Ref. [30], this scenario might not 
be compatible with semi-local duality, when the number of colors is sufficiently high. Besides, we should note that, in sharp 
contrast with the work here, it is not clear whether the Ttn amplitudes used in Ref. [70] contain or not all the leading 1/Nc 
contributions and thus, the analysis of the behavior of the resonance properties when Nc is larger than 3 is not fully meaningful. 

The qualitative Nc behavior of the two resonances found in the scalar-isoscalar sector is substantially different when one 
looks into the results of the fits B.2 and B.1-2 (middle panels of Figs. |9] and [TOli. There, the role played by the Nc — 3 G and 
scl resonances is interchanged. Thus, the qq component of the a seems to be absent, and the /o(600) pole becomes wider and 
moves into the third quadrant above Nc > 5. Indeed, now the behavior showed by this pole is quite similar to that displayed for 
the (7 in the right panel of Fig. 4 of Ref. 1701 . On the other hand, as can be appreciated in the middle panels of Fig. [TOl the second 
I = J = resonance now becomes the scalar SRA pole included in our amplitudes. Presumably, at high Nc it would become 
more delta-function-like, as the p pole would, and it would likely provide the needed cancellations with the contribution of this 
latter resonance in the elastic n'^n^ amplitude OOll . 



A recent study II72I1 describes an accidental symmetry of the full Regge tower of radially excited 0++ states, M — an + m 



a' 



Remarkably, the states generating doublets with excited pion states are /o(600) ^ ;ro(140), /o(1370) O ;ro(1300), /o(1710) ^ 
7Cb(1800), /o(2100) O TZa[2QlQ) and /o(2330) O ;^o(2360), whereas the other scalar states /o(980), /o(1500), /o(2020) 
and /o(2200) are not degenerate with other mesons with light u and d quarks. Arguments have been put forward as to 
identify the /o(980) as a would-be glueball in the large-A^c limit. As it is well known, glueballs are more weakly cou- 
pled to mesons, (?{\/Nc), than other mesons, ff{\/^/Nc). This is supported by the rather small width ratio which yields 
r//r<j -- {g}j,Tim^f) / [ginK^n],) ^ \/Nc, and for nic ^ 0.8 MeV a ratio gcTm/gfim ^ VW is obtained. 

Our analysis can be improved along several lines. Firstly, one might extend the chiral analysis to include two-loop results. 
Secondly, our conclusions might be modified when coupled channels incorporating KK effects are taken into account. We have 
given arguments why we do not expect that this might be important as long as we remain in the sub-threshold region, where all 
KK effects should be included as l/Nc corrections to the counter- terms. 

VII. CONCLUDING REMARKS 

We summarize the conclusions of this work. First, we have constructed tzk amplitudes that fulfill exact elastic unitarity, 
account for one-loop ChPT contributions and include all leading terms, within the SRA, in the 1 /Nc expansion. These amplitudes 
have been successfully fitted to 7 = 7 = 0, 7 = 2, 7 = and 7 = 7=1 phase shifts. Next, we have looked for poles in the SRS of 
the amplitudes, and discussed their properties. Since all leading 1 /Nc terms are taken into account, this scheme is much more 
appropriated to discuss the Nr dependence of the (7 and p masses and widths than previous ones, where the leading 1 /Nc terms 
appearing beyond ff{p'^) [la, |23[ or 6'{p^) \1V\ were neglected. The recent work of Ref. [^29] does not identify correctly the 
leading 1 /Nc term, and hence the conclusions of this reference in the large-A/c limit need some revision. 



■^' Note however, that such a statement requires an extrapolation of the Regge behavior to somehow low energies. For that end, the authors of Ref. \iO\ use 
Regge trajectories for the variable ( V — 2m^ — f/2) to ensure that the imaginary part of the extrapolated Regge amplitude vanishes at threshold. 
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Robust conclusions are drawn in the case to the p resonance, and we confirm here that it is a stable meson in the large-A^c 
limit, as pointed out by other authors in the past. In the scalar-isoscalar sector, the overall scenario looks like somehow less 
predictive, since we cannot firmly conclude whether or not the Nc — 'ifo (600) resonance completely disappears at large Nc or it 
has a sub-dominant component in its structure, which would become dominant when the number of colors gets sufficiently high. 
Unfortunately, this depends on the chosen procedure (A, B.l and B.3 or B.2 and B.1-2) to fit the phase-shift data. However, it 
becomes clear the predominant di-pion component of this pole for Nc = 3, and that the SRA delta-function-like pole, that always 
appears in the Nc ^ I limit, might help to keep the whole scheme compatible with semi-local duality. Nevertheless, this needs 
to be quantitatively checked elsewhere. 

Acknowledgments 

We thank J.R. Pelaez, J. Ruiz de Elvira and J.J. Sanz-Cillero for useful communications. This research was supported by 
DGI and FEDER funds, under contracts FIS2008-01143/FIS, FRA.2007-60323 and the Spanish ConsoHder-Ingenio 2010 Pro- 
gramme CPAN (CSD2007-00042), by Junta de Andalucia contract FQM0225, by Generalitat Valenciana under contracts PROM- 
ETEO/2008/069 and PROMETEO/2009/0090, and it is part of the European Community-Research Infrastructure Integrating 
Activity "Study of Strongly Interacting Matter" (acronym HadronPhysics2, Grant Agreement n. 22743 1), under the Seventh EU 
Framework Programme. A. P. acknowledges the support of the Alexander von Humboldt Foundation. 



[1 
[2 

[3 

[4; 

[5 
[6 

[7 

[8: 

[9 
[10 
[11 
[12; 
[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 

[22: 

[23 
[24 
[25 
[26 
[27: 
[28 
[29: 
[30 
[31 
[32 
[33 
[34: 
[35 
[36 

[37: 

[38 
[39: 
[40 



M. H. Johnson and E. Teller, Phys. Rev. 98, 783 (1955). 

I. Caprini, G. Colangelo, and H. Leutwyler, Phys. Rev. Lett. 96, 132001 (2006), hep-ph/05 12364. 

R. Kaminski, J. R. Pelaez, andP. J. Yndurain, Phys. Rev. D77, 054015 (2008), 0710.1150 [hep-ph]. 

R. Garcia-Martin, R. Kaminski, J. Pelaez, and J. de Elvira (201 1), 1 107.1635 [hep-ph]. 

E. Klempt and A. Zaitsev, Phys. Rept. 454, 1 (2007), 0708.4016 [hep-ph]. 

J. R. Pelaez, Phys. Rev. Lett. 92, 102001 (2004), hep-ph/0309292. 

G. 't Hooft, Nucl. Phys. B72, 461 (1974). 

E. Witten, Nucl. Phys. B160, 57 (1979). 

G. Ecker, J. Gasser, A. Pich, and E. de Rafael, Nucl. Phys. B321, 311 (1989). 

G. Ecker, J. Gasser, H. Leutwyler, A. Pich, and E. de Rafael, Phys. Lett. B223, 425 (1989). 

S. Weinberg, Physica A96, 327 (1979). 

J. Gasser and H. Leutwyler, Ann. Phys. 158, 142 (1984). 

J. Gasser and H. Leutwyler, Nucl. Phys. B250, 465 (1985). 

V. Cirigliano et al., Nucl. Phys. B753, 139 (2006), hep-ph/0603205. 

J. Bijnens, E. Gamiz, E. Lipartia, and J. Prades, JHEP 04, 055 (2003), hep-ph/0304222. 

E Sannino and J. Schechter, Phys. Rev D52, 96 (1995), hep-ph/9501417. 

M. Harada, E Sannino, and J. Schechter, Phys. Rev. D54, 1991 (1996), hep-ph/95 11335. 

A. Pich (2002), hep-ph/0205030. 

J. Nieves and E. Ruiz Arriola, Phys. Lett. B679, 449 (2009), 0904.4590 [hep-ph]. 

Z. X. Sun, L. Y. Xiao, Z. Xiao, and H. Q. Zheng, Mod. Phys. Lett. A22, 711 (2007), hep-ph/0503195. 

J. R. Pelaez and G. Rios, Phys. Rev. Lett. 97, 242002 (2006), hep-ph/06 10397. 

J. Nieves and E. R. Arriola, Phys. Rev. D80, 045023 (2009), 0904.4344 [hep-ph]. 

P Di Vecchia and G. Veneziano, Nucl.Phys. B171, 253 (1980). 

E. Witten, Annals Phys. 128, 363 (1980). 

C. Rosenzweig, J. Schechter, and C. Trahem, Phys.Rev. D21, 3388 (1980). 
R. Kaiser and H. Leutwyler, Eur Phys. J. C17, 623 (2000), hep-ph/0007101. 
N. Beisert and B. Borasoy, Phys. Rev. D67, 074007 (2003), hep-ph/0302062. 

M. Albaladejo, J. A. Oiler, and L. Roca, Phys. Rev D82, 094019 (2010), 1011.1434 [hep-ph]. 

Z.-H. Guo and J. A. Oiler (201 1), 1 104.2849 [hep-ph]. 

J. Ruiz de Elvira, J. R. Pelaez, M. R. Pennington, and D. J. Wilson (2010), 1009.6204 [hep-ph]. 

B. Martin, D. Morgan, and G. Shaw, Pion-pion interactions in particle physics, vol. 101 (Academic press London, 1976). 
E J. Yndurain (2002), hep-ph/02 12282. 

J. Nieves and E. Ruiz Arriola, Nucl. Phys. A679, 57 (2000), hep-ph/9907469. 

R. Garcia-Martin, R. Kaminski, J. Pelaez, J. Ruiz de Elvira, and F. Yndurain, Phys.Rev. D83, 074004 (2011), 1102.2183 [hep-ph]. 

J. R. Pelaez and E J. Yndurain, Phys. Rev. D68, 074005 (2003), hep-ph/0304067. 

J. R. Pelaez and E J. Yndurain, Phys. Rev. D69, 114001 (2004), hep-ph/0312187. 

S. L. Adler and E J. Yndurain, Phys. Rev. D75, 1 16002 (2007), 0704. 1201 [hep-ph]. 

D. Black, A. H. Eariborz, E Sannino, and J. Schechter, Phys. Rev. D59, 074026 (1999), hep-ph/9808415. 
V. Cirigliano, G. Ecker, H. Neufeld, and A. Pich, JHEP 06, 012 (2003), hep-ph/030531 1. 

E Giacosa, T Gutsche, V. E. Lyubovitskij, and A. Eaessler, Phys. Rev. D72, 094006 (2005), hep-ph/0509247. 



27 

[41] Z.-H. Guo and J. J. Sanz-Cillero, Phys. Rev. D79, 096006 (2009), 0903.0782 [hep-ph]. 

[42] Z. H. Guo, J. J. Sanz Cillero, and H. Q. Zheng, JHEP 06, 030 (2007), hep-ph/0701232. 

[43] A. Pich, I. Resell, and J. J. Sanz-Cillero, JHEP 02, 109 (201 1), 101 1.5771 [hep-ph]. 

[44] Z.-H. Guo and P Roig, Phys. Rev. D82, 113016 (2010), 1009.2542 [hep-ph]. 

[45] G. Ecker, Acta Phys. Polon. B38, 2753 (2007), hep-ph/0702263. 

[46] R J. Oilman and H. Harari, Phys. Rev. 165, 1803 (1968). 

[47] S. Weinberg, Phys. Rev. 177, 2604 (1969). 

[48] S. Weinberg, Phys. Rev. Lett. 65, 1177 (1990). 

[49] G. Ecker and C. Zauner, Eur. Phys. J. C52, 315 (2007), 0705.0624 [hep-ph]. 

[50] D. Toublan, Phys. Rev. D53, 6602 (1996), hep-ph/9509217. 

[51] K. N. et al., J. Phys. G37, 075021 (2010). 

[52] C. Garcia-Recio, L. Geng, J. Nieves, and L. Salcedo, Phys.Rev. D83, 016007 (2011), 1005.0956 [hep-ph]. 

[53] E. Megias, E. Ruiz Arriola, L. L. Salcedo, and W. Broniowski, Phys. Rev. D70, 034031 (2004), hep-ph/0403139. 

[54] J. Nieves and E. Ruiz Arriola, Phys. Lett. B455, 30 (1999), nucl-th/9807035. 

[55] I. Rosell, J. J. Sanz-Cillero, and A. Pich, JHEP 0408, 042 (2004), hep-ph/0407240. 

[56] I. Rosell, J. J. Sanz-Cillero, and A. Pich, JHEP 0701, 039 (2007), hep-ph/06 10290. 

[57] A. Pich, I. Rosell, and J. J. Sanz-Cillero, JHEP 0807, 014 (2008), 0803.1567 [hep-ph]. 

[58] I. Rosell, P Ruiz-Femenia, and J. Portoles, JHEP 0512, 020 (2005), hep-ph/05 10041. 

[59] J. Portoles, I. Rosell, and P Ruiz-Femenia, Phys.Rev. D75, 114011 (2007), hep-ph/06 11375. 

[60] J. Nieves and E. Ruiz Arriola, Phys. Rev. D64, 116008 (2001), hep-ph/0104307. 

[61] B. Ananthanarayan, G. Colangelo, J. Gasser, and H. Leutwyler, Phys. Rept. 353, 207 (2001), hep-ph/0005297. 

[62] F J. Yndurain, R. Garcia-Martin, and J. R. Pelaez, Phys. Rev. D76, 074034 (2007), hep-ph/0701025. 

[63] P Estabrooks and A. D. Martin, Nucl. Phys. B79, 301 (1974). 

[64] S. D. Protopopescu et al., Phys. Rev. D7, 1279 (1973). 

[65] W. Hoogland et al., Nucl. Phys. B126, 109 (1977). 

[66] M. Losty, V. Chaloupka, A. Ferrando, L. Montanet, E. Paul, D. Yaffe, J. Zieminski, et al.. Nuclear Physics B 69, 185 (1974), ISSN 

0550-3213. 
[67] G. Amoros, J. Bijnens, and P Talavera, Nucl. Phys. B585, 293 (2000), hep-ph/0003258. 
[68] J. Bernabeu and J. Prades, Phys. Rev. Lett. 100, 241804 (2008), 0802.1830 [hep-ph]. 
[69] H. Leutwyler, AIP Conf. Proc. 1030, 46 (2008), 0804.3182 [hep-ph]. 
[70] Z.-Y. Zhou and Z. Xiao, Phys. Rev. D83, 014010 (2011), 1007.2072 [hep-ph]. 
[71] N. A. Tornqvist, Z. Phys. C68, 647 (1995), hep-ph/9504372. 
[72] E. Ruiz Arriola and W. Broniowski, Phys. Rev. D81, 054009 (2010), 1001.1636 [hep-ph]. 



